Historic,  archived  document 

Do  not  assume  content  reflects  current 
scientific  knowledge,  policies,  or  practices 


XEROX   COPY     DO  NOT  REMOVE 


A99.9  F764Un 


UNITED   STATES    DEPARTMENT  OF  AGRICULTURE 
FOREST  SERVICE 

mmuomwa  forest  &  uhe  eiperimem  station 

OCDEN  UTAH 


USDA  Forest  Service 

Research  Note  INT-98  1969 

COMPUTING  PROCEDURE  FOR  GRAND  FIR  SITE 
EVALUATION  AND  PRODUCTIVITY  ESTIMATION 

Albert  R.  Stage1 


ABSTRACT 


Two  subroutine  programs  to  calculate  productivity  of 
grand  fir  sites  from  measurements  of  tree  height,  crown 
ratio,  and  age  are  described.    A  method  for  obtaining  a 
graphical  estimate  of  site  quality  is  also  provided  when 
machine  calculation  is  not  appropriate.    These  procedures 
evaluate  site  on  scales  directly  proportional  to  the 
potential  of  the  site  for  growing  grand  fir.  This 
potential  is  measured  in  units  of  cubic  feet  of  the  net 
current  annual  increment  of  the  entire  stand  at  the  staae 
of  development  at  which  mean  annual  increment  culminates. 

A  nonlinear,  implicit  expression  for  relative  productivity  (Q)  of  grand  fir  sites 
was  introduced  in  a  previous  paper  (Stage  1966)  along  with  a  procedure  for  rating  pro- 
ductivity in  cubic  feet  per  acre  per  year.     When  relative  productivity  of  the  site  is 
estimated  from  data  for  individual  trees,  there  is  considerable  variation  of  these 
estimates  due  to  the  relative  dominance  of  the  trees  within  the  stand.    Much  of  this 
variation  can  be  removed  by  adjustments  based  on  the  ratio  of  live  crown  to  total 
height  of  the  tree.     The  purpose  of  this  Note  is  to  provide  computing  algorithms  to 
help  forest  managers  evaluate  site  and  productivity  using  these  procedures. 

The  expression  for  relative  productivity  is 

Q  =  bsl  (h*  -  4.5)b2  (1  -  e"b3Q(t+b5))bl+  CD 

in  which 


t  =  age  at  breast  height  in  years 

c  =  live  crown  length  in  percent  of  total  height 

h*=  total  height  in  feet  (h)  adjusted  for  the  departure  of  crown  length 
from  the  standard  cro\vn  percentage  for  a  tree  of  this  age. 
[In  (h*  -  4.5)  =  In  (h  -  4.5)  -  18.786  (In  c  -  3.747474  -  8 . 218305/t) /t] 


Principal  Mensurationist ,  stationed  in  Moscow,  Idaho,  at  Forestry  Sciences  Lab- 
oratory, maintained  in  cooperation  with  University  of  Idaho. 


bx=  +3.368X10-5 
b2=  +2.00 
b3=  +0.018167 


t>4=  -4.00 
b5=  -2.70 


A  listing  of  a  FORTRAN  subroutine  to  solve  expression  (1)  is  given  in  the  Appendix. 
This  program  is  derived  from  a  general  purpose  program  to  solve  general  nonlinear 
equations,  such  as  x  =  f(x)   (Lance  1960;  Wegstein  1960)  and  is  executed  by  the 
statement 

CALL  RPGF  (Q,H,T,C,IE) 

in  which  the  arguments  are  variables  defined  as  follows: 

Q  =  relative  productivity  as  computed  by  RPGF 
H  =  height  in  feet 

T  =  age  in  rings  at  4.5  feet  above  ground 

C  =  live  crown  length  in  percent  of  total  height 

IE=  an  integer  variable  that  should  be  zero  on  return  from  RPGF.  Nonzero 
values  of  IE  indicate  one  of  several  possible  error  conditions  detected 
by  RPGF.     These  error  codes  are  listed  in  the  Appendix. 

Q  can  also  be  estimated  from  table  1  and  figure  1.     First,  the  total  tree  height 
is  adjusted  for  crown  ratio  by  the  equation 

h*  =  (h  -  4.5)   (multiplier  from  table  1)  +  4.5. 

Then,  Q  can  be  read  at  the  point  corresponding  to  h*  and  age  at  breast  height  in  figure 
1  just  as  you  use  conventional  site  index  curves. 

Table  1. --Height  multiplier  for  crown  ratio  adjustment 
in  grand  fir  site  evaluation 
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Figure  1. — Curves  to  estimate  Q  from  height  (adjusted  for  crown  ratio)  and  age  at 

breast  height. 

To  convert  Q,  the  estimate  of  relative  productivity,  to  productivity  in  net  cubic 
feet  per  acre  per  year  requires  an  additional  subprogram,  the  FORTRAN  function  given 
below: 


FUNCTI  C.N   CUPGF  {  S  ) 

8   =   89.761/S  -  22.701 

A   =  B   *(0.5   +   SQRT(C25   +   12./(S*B)  )  I 
CUPGF   =   22228. *EXP(-B/A) /{  A  +    12. /S) 
RETURN 
END 


This  subprogram  is  written  as  a  function  to  be  executed  by  a  statement  such  as 

PROD  =  CUPGF(Q) 

in  which  the  variable  PROD  will  be  assigned  the  value  of  the  productivity  computed  by 
the  function  CUPGF  using  Q  (from  RPGF)  for  its  argument. 


The  following  tabulation  shows  estimates  of  productivity  associated  with  values 
of  Q  from  0.6  to  1.5. 

Q  value  Productivity 

(Cu. ft . /acre/year) 

0.6  56 

.7  67 

.8  79 

.9  91 

1.0  105 

1.1  118 

1.2  133 

1.3  149 

1.4  166 

1.5  184 


REFERENCES 

Lance,  G.  N. 

1960.     Numerical  methods  for  high-speed  computers.     Iliffe,  London.     166  pp. 
Stage,  A.  R. 

1966.     Simultaneous  derivation  of  site-curve  and  productivity  rating  procedures. 
Soc.  Amer.   Forest.  Proc.   1966:  134-136. 

Wegstein,  J. 

1960.     Algorithm  2.     Commun.  Assoc.  Comp.  Mach.  3(2):  74. 


4 


APPENDIX 


Error  Messages  (IE) 

The  program  RPGF  will  find  a  solution  for  most  realistic  combinations  of  height 
and  age  of  dominant  trees.     However,  the  algorithm  may  fail  as  a  result  of  errors  in 
the  tree  data  or  poor  selection  of  trees  for  site  evaluation.     If  so,  the  variable  IE 
will  be  nonzero  on  return  from  RPGF. 

If  IE  is  positive,  then  the  20  programed  iterations  were  not  sufficient  to  attain 
the  level  of  precision  specified  in  the  accuracy  test  for  termination.     The  solution 
may  not  be  correct,  even  when  less  stringent  accuracy  criteria  are  applied. 

If  IE  =  -2,  then  the  derivative  of  the  growth  function  equaled  zero  for  some  trial 
estimate  of  Q.     This  condition  could  result  from  the  starting  values  built  into  RPGF, 
or  from  accumulation  of  round-off  errors  in  the  machine  arithmetic.     (In  the  author's 
experience  with  data  from  several  thousand  trees,  this  error  condition  has  never 
occurred . ) 

If  IE  =  -3,  several  data-based  errors  may  be  at  fault.     A  value  for  age  of  less 
than  one  is  the  most  likely  cause.     Otherwise,  the  error  condition  might  be  due  to 
unreasonable  combinations  of  age  and  height  that  cause  divergence  of  the  algorithm. 

Information  for  Programers 

These  programs  use  the  following  library  routines: 

ALOG(X)  =  natural  logarithm  of  X 

DLOG(X)  =  double  precision  version  of  ALOG(X) 

EXP(X)  =  exponential  function  of  X 

DEXP(X)  =  double  precision  version  of  EXP(X) 

DABS(X)  =  absolute  value  of  X,  double  precision. 

Computers  with  floating  point  words  longer  than  32  bits  would  not  require  the 
double  precision  arithmetic  that  is  used  in  the  programs  shown  on  page  3  and  below. 
To  convert  to  single  precision,  the  names  of  the  functions  must  be  changed  appro- 
priately, the  DOUBLE  PRECISION  statement  removed,  and  the  "D"  changed  to  "E"  in  the 
exponent  of  the  constants  (e.g.,  instructions  91,  991,   1  and  3,  below). 

FORTRAN  subroutine  to  calculate  relative  productivity   

  ,    k  — - 

SUBROUTINE   RPGF(    Q»HjTtCRi  IGR  ) 

jOU3L£   P'^iC  I  SiKN   A  ,  ii ,  C  ,  D  ,  X  ,  BT  ,  OET  ,  TOL  ,  VAL,  BRAG,  ALGCR 
DIMENSION  C(d) 

0  AT  A   C  /-1C. 29862,    2.0,    0.016167,    -2.7    ,    2.0  / 
DIV   =  0.5 

r  t.R  =  o 
icr  =  i 

COR   =  0. 

1  =  c 

.U   =-C(3)*(T        C<4)  ) 
IF    (    T  -   90.    )    91,91 ,92 
TOL   =   0. GOO 

GD  TC  9  5 

T;jL   =   ALCIG(    0.3*H/T  ) 

ftLGCW    =   AL0G(H-4.«3)    -    I  ALGG  (CR  )-3. 747474-8.  218305/T ) *18. 786/T 

(con.  next  vage) 


FORTRAN  subroutine  to  calculate  relative  productivity  (con.) 
TS  =  TOL 

99    IF    (    TOL    .GT.    20.)    GO   TO  101 

UET   =   BT*UEXP(    TGL  ) 
IF    (    DET    )  991,100,100 

991  BR AC   =    1.000  -  QEXP(DET) 
IF    ( BR AC )  100,100,992 

992  TGL  =  C(l)    +  C(2)-(ALGCR  -  C ( 5 ) *DL0G ( BR AC )  ) 
1   =    1   +  1 

996   GG  TG    {    21,22,23) , I GT 

21  X   =  TOL 

A   =   X  -  TS 
8  =  -A 
I  GT   =  2 
•   GG   TC  99 

22  VAL   =   X  -  TGL 

C  START    ITERATION  LOOP 

IF    (    7 AL )  1,7,1 
C  EQUATION    IS  NOT   SATISFIED   BY  X 

1  B   =   b  /VAL  -  l.OOC 
IF(B)    2, «, 2 

C  ITERATION    IS  POSSIBLE 

2  A   =  A/3 

X    =    X    +  A 
L>    =  VAL 
TGL   =  X 
IoT  =  3 
GG  TG  99 

23  VAL   =   X  -  TOL 

C  TEST  ON   SATISFACTORY  ACCURACY 

0  =  DABS(X) 
TOL   =  0.0008*0 

4  IF    (    UABS(A)-rGL)  5,5,6 

5  I F ( DABS ( VAL )- 1 • 0D0*TCL )  7,7,6 

6  IF    {    1-20    )  1,1,10 

C  ENU  OF    ITERATION  LOOP 

7  3   =  DEXP(X) 
RETURN 

C  ERROR   RETURN   IN  CASE  OF   ZERO  DIVISOR 

8  I ER   =  -2 
RETUR  1 

C  NL  CONVERGENCE   AFTER   20    ITERATION  LOOPS.    ERROR  RETURN 

10    I £R   =  I 
RETURN 

C  ERROR   CONDITION  ON    INITIAL  ITERATION 

12    I  ER   =  -3 

RETURN 
14   I  E9   =  -4 

RETURN 

100  GU  TG    ( 1 ?  ,  12  ,  102 )  , I GT 

101  GO  TG    (    14, 14, 102 ) , IGT 

102  X   =   X  -  A 

A    =   D  I  V  *  A  /  P 
D I  V   =   D I  V  *  D  I  V 

1  =   1   +  1 
GO   TO  2 

E  ID 
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